Landau level spectroscopy of ultrathin graphite layers 
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Far infrared transmission experiments are performed on ultrathin epitaxial graphite samples in a 
magnetic field. The observed cyclotron resonance-like and electron-positron-like transitions are in 
excellent agreement with the expectations of a single-particle model of Dirac fermions in graphene, 
with an effective velocity of c = f.03xf0 6 m/s. 

PACS numbers: 71.70.Di 76.40.+b 78.30.-j 78.67.-n 



i 

-a 
c 

o 
o 



(N 
> 

CO 

in 
o 
\o 
o 

-I— > 
c3 



i 

c 

o 
o 



X 



The electronic properties of graphite have recently be- 
come the center of considerable attention, following ex- 
periments on graphite monolayers (graphene) [lj and epi- 
taxial graphene [2|, which led to the discovery of an un- 
usual sequence of quantum Hall effect states 0, H[ and 
an energy-dependent mass. The considerable interest 
in two-dimensional graphite is fuelled by its particular 
band structure and ensuing dispersion relation for elec- 
trons, leading to numerous differences with respect to 
"conventional" two-dimensional electron systems (2DES) 
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0, EH, on single crystal (4H) SiC. These epitaxial 
graphene structures are routinely characterized using low 
energy electron diffraction, Auger electron spectroscopy, 
X-ray diffraction, scanning tunnelling microscopy and 
atomic force microscopy. The results of these measure- 
ments in combination with angular resolved photoclcc- 
tron spectroscopy and transport data indicate that the 
graphitized part of this type of structure consists of a few 
(3-5) graphene layers 0, EH- We investigated two such 
(unpatterned) structures, with dimensions of about 4 x 



111 I X 21 ] - The band structure of 4 mm , which both show a similar behavior. 



graphene is considered to be composed of cones located at 
two inequivalent Brillouin zone corners at which the con- 
duction and valence bands merge. In the vicinity of these 
points the electron energy depends linearly on its mo- 
mentum: E(~p) — ±c| p |, which implies that free charge 
carriers in graphene are governed not by Schrodinger's 
equation, but rather by Dirac's equation for zero rest 
mass particles, with an effective velocity c, which replaces 
the speed of light. With the application of an external 
magnetic field, the Dirac energy spectrum evolves into 
Landau levels with energies given by 



E n = sgn(n)c\/2ehB\n\ — sgn(n)Ei 



(1) 



where n scans all positive (for electrons) and negative 
(for holes) integers and - very importantly - zero. E\ 
may be understood as a characteristic energy introduced 
by the magnetic field. The square root dependence on B 
and Landau level index n is in stark contrast to "conven- 
tional" 2D electrons, where E n = (n+^)heB /m, (n > 0), 
and the Landau levels are equally spaced. 

The unusual sequence of quantum Hall effect states 
and an energy-dependent electron effective mass [H, 0], 
found in magneto-resistance measurements, are consis- 
tent with the model of Dirac particles. Here we report a 
magneto-spectroscopy study of the optical properties of 
ultrathin epitaxial graphite layers, in which we directly 
probe the dependence of the energy of electrons on their 
momentum. 

The experiments were performed on graphene layers 
grown in vacuum by the thermal decomposition method 



The far infra-red transmission of the samples was mea- 
sured, at a temperature of 1.9 K, as a function of the 
magnetic field B. A Si bolometer was placed directly be- 
neath the sample to detect the transmitted radiation. 
The light (provided and analyzed by a Fourier transform 
spectrometer) was delivered to the sample by means of 
light-pipe optics. All experiments were performed with 
non-polarized light, in the Faraday geometry with the 
wave vector of the incoming light parallel to the mag- 
netic field direction. The transmission spectra were nor- 
malized by the transmission of the substrate and by the 
zero-field transmission, thus correcting for magnetic field 
induced variations in the response of the bolometer. The 
SiC substrate used was completely opaque for energies 
between 85 meV and about 200 meV, which limited the 
range of our investigation. 

The main experimental finding consists of several ab- 
sorption lines visible in the spectra. A representative 
transmission spectrum for sample 1, at 0.4 T, is shown 
in Fig. [1] These lines evolve spectacularly with the mag- 
netic field. Two main lines are shown in Fig. [2] for sev- 
eral values of the field. As shown in Fig. [3l their en- 
ergies, plotted as a function of the square root of the 
magnetic field, trace perfect straight lines, in excellent 
agreement with eq. [T] Also shown in this figure are the 
energy positions of the two other lines. Experiments per- 
formed in a tilted configuration show that the position 
of the transition line (filled symbols in Fig. [3]) depends 
only on the component of the magnetic field perpendicu- 
lar to the sample plane. A complicated structure at still 
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FIG. 1: (Color online) Relative transmission trace at 0.4 T 
and 1.9 K shows 4 distinct transitions. The assignations are 
(see text) A: Li — > L2, B: Lo — » Li(L_i — > Lo), C : L_2 — > 
Li(L_i — > La), D: L_3 — > Li(L-2 — * L3). The inset shows 
a schematic of the evolution of Landau levels with applied 
magnetic field, and possible optical transitions. 



lower energies, not shown in the figure, moving slowly to 
higher energies with magnetic field, was also observed. 
Unusually, the intensities of the two main lines increase 
markedly with increasing magnetic field, with the lower- 
energy line always remaining stronger. 

Cyclotron resonance in graphite has been studied ex- 
perimentally [lH and theoretically 14, 15 1. These ex- 
periments showed a linear dependence of the cyclotron 
frequency on the magnetic field, with an effective mass 
of 0.058nio. Our results are best described using the pre- 
dictions of a simple single-particle (Dirac) model for a 
graphene layer, and we will use this language in the fol- 
lowing paragraphs. To facilitate discussion, we sketch the 
graphene Landau levels and possible transitions between 
them in the inset to Fig. Q] 

Thus, we assign the strongest line to the transitions 
to and from the lowest Landau level. Note that, since 
the conduction and valence band states in graphene are 
built from the same atomic orbitals, the positive and 
negative branches of the dispersion relation are identi- 
cal. The ensuing symmetry means that the Lq — ► L\ 
and L_i — > Lq transitions are indistinguishable in an 
experiment using unpolarized radiation. A straight line 
fit of the points corresponding to this transition using 
the expression E = E\ = cy2ehB yields a very accurate 
value for c, the velocity of electrons in graphene. This 



is found for both samples to be (1.03 ±0.01) xl0 6 m/s, 
consistently with transport measurements [H, 0] • 

The slopes of the other lines traced in Fig. [3j starting 
from the highest energy transition, scale exactly as (V3+ 
\/2) : (V2+ 1) : 1 : [y/2 - 1), allowing these lines to be 
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FIG. 2: (Color online) Development of two main transmis- 
sion lines with the magnetic field. Upper panel shows the line 
marked C in Fig[T] corresponding to L_i — » La(L_2 — * Li) 
transitions, lower panel shows the line marked B, correspond- 
ing to Lq — > Li(L_i — + Lq) transitions: solid lines are for 
sample 1, dashed lines for sample 2. 



assigned to transitions L_2 — > £3(^-3 — * L2), L_i — > 
L 2 (L_ 2 -> Li), L a -> Li(L_ x -> Lq), and L x -> L 2 , 
respectively, as shown in the figure. 

The fact that transitions involving the Lq Landau level 
are visible at such low magnetic fields places an upper 
limit on the electron concentration in the observed layer. 
The observation of the Lq — > L\(L-\ — > Lq) line implies 
the existence of unpopulated states at least on the L\ 
level. This line is clearly observed at fields B w 0.15 T, 
and therefore the L\ level can be fully populated only 
when B < 0.15 T. Thus n < 2.1 x 10 10 cm" 2 (where we 
take into account the 2- and 4-fold degeneracy of the Lq 
and L\ electronic Landau levels, respectively). This is 
also consistent with the disappearance of line A (L\ — > 
L2) when the Li level is depopulated by the magnetic 
field (see Fig.gJ). 

We now turn our attention to the strength of the tran- 
sitions. As may be seen in Fig. [2l both the main transi- 
tions gain in intensity with increasing magnetic field. To 
better visualize this trend, we plot the integrated inten- 
sity(area under the dip in the relative transmission) for 
sample 1, as a function of the square root of the mag- 
netic field (Fig. 2]). Sample 2 shows the same behavior 
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FIG. 3: The observed transitions, together with their assign- 
ments, plotted versus B±. The filled symbols are data ob- 
tained for the sample tilted with respect to the direction of B 
by an angle of 50°. The dashed lines are all calculated with 
the same characteristic velocity £ = 1.03 x 10 6 m/s. 




FIG. 4: Area of the main peak observed in the experiments 
- the Lq — > L\(L-\ — > Lo) transition - plotted against the 
square root of the magnetic field. The line is traced using the 
expression 



but smaller values of the intensity. 

The relative transmission of a sheet of conducting elec- 
trons between vacuum and a dispersionless polar medium 
with a refractive index k, for unpolarized radiation and 
in the limit of weak absorption, may be written as (see 



e.g. 



IE 03) 



T(w,B)«l-/3 



Re(a xx (uj,B)) 



e c 



where cr xx (u>, B) is a diagonal element of the optical con- 
ductivity tensor, eo is the vacuum permittivity, c is the 
speed of light in vacuum and (3 = (n 2 + 3)/2(k 2 + 1) = 
0.63 for SiC, where k—2.6 (from the substrate transmis- 
sion) . The optical conductivity of the 2D electrons may 
be written using the Kubo formalism [l8| and taking into 
account the prop erties of the graphene Landau level wave 
functions 
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where E m>n are the transition energies between levels m 
and u,Gb = eB/h is the Landau level degeneracy, f m , f n 
are the occupancies of the relevant Landau levels and 
the selection rules for the optically active transitions are 
given by M m , n = (c 2 /p)S\ m \ ! \ n \±i, with p = 2 for m or n 
=0 and 4 otherwise. The summation is performed over 
all Landau levels m, n, and the fourfold degeneracy of 
each Landau level has already been accounted for. 

The integrated transmission for a single transition be- 
tween a completely filled (Lo) & n( A a completely empty 
(Li) Landau level, using the above expression (for 



= /3 



e 2 c 



linewidths 7 <C E m ^ n ) may be written as: 

1 f e 3 c 2 B 

1(B) = — / Re{o xx {u))du w — , . 

e c J £qcEi 2e hc 

where E\ is the characteristic energy introduced earlier. 

The above equation gives a rough estimate of the in- 
tensity of the strongest transition, in the range of high 
magnetic fields where the Fermi energy is pinned to the 
Lq level. This is due to the fact that the decreasing inten- 
sity of the Lq — > L\ transition is compensated by the cor- 
responding increase of the strength of the superimposed 
L_i — > Lo transition. At low magnetic fields, where the 
L\ level is not completely empty, the observed oscilla- 
tor strength decreases, disappearing when the L\ is fully 
populated. In spite of the rather crude approximation, 
Fig. 0] indeed shows that the observed transition follows 
the expected trend; the good agreement of the absolute 
measured and calculated values is another factor sup- 
porting the picture of a single, possibly inhomogeneous, 
graphene layer (see discussion in following paragraphs). 

Several notable differences emerge between Dirac elec- 
trons and conventional two-dimensional electron systems 
(2DES). As we have shown, transitions between adjacent 
Landau levels in graphene occur at markedly different en- 
ergies (for example Lq — > L\ and L\ — > L2, Figs.Q]and[3]). 
For a standard 2DES, transitions between such pairs of 
Landau levels all have the same energy. More striking, 
a different class of transitions, with no counterpart in 
a standard 2DES, is observed in graphene and involves 
those from hole (n < 0) to electron (n > 0) states (e.g. 
our L-i — » L2 and L_2 — * L\ transitions). These are 
the particle-antiparticle creation and annihilation events 
in the Dirac formalism. 
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Since some of the observed transitions are analogues 
of cyclotron resonance, it is tempting to look at them 
in a semi-classical context, using the concept of an effec- 
tive mass. While for a 2DES with a quadratic dispersion 
law there is a coincidence between classical and quantum 
mechanical solutions of the optically active response in a 
magnetic field, this does not hold for graphene. The clas- 
sically derived cyclotron excitation Ec in this system is 
Ec = heB/(E/c 2 ) [H, where E is the electron energy and 
(E/c 2 ) stands for the electron mass. Although the effec- 
tive rest mass of the electrons in graphene is zero, their 
energy- and magnetic field-dependent cyclotron mass can 
be followed down to the lowest energies (~ 7 meV in our 
case), giving a lowest observed value of 0.0012 mo- 

Having demonstrated the presence of zero effective rest 
mass Dirac fermions in the investigated structure, let 
us now consider the following points: (i) linear disper- 
sion is characteristic of a single graphene layer, while a 
graphene bilayer [2(| 21 1 is found to exhibit parabolic 



dispersion; (ii) transport measurements performed on a 
mesoscopic sample patterned on the same wafer as our 
selected sample, which also show the unusual Berry's 
phase of n observed in graphene, give a concentration 
of ~ 4 x 10 12 cm -2 , while our results point to a concen- 
tration two orders of magnitude smaller. 

It is believed that electric transport is dominated by 
the interface layer, which has a high electron concen- 
tration due to the built-in electric field caused by the 
surface charge 0, [Hj]. Transmission measurements, on 
the other hand, probe the whole sequence of layers, in- 
cluding those further away from the interface which have 
lower electron concentrations. It is possible that the ob- 
served Dirac spectrum originates from a single graphene 
layer "floating" above a SiC substrate covered with other 
graphitic layers [23] . The previously mentioned very low- 
energy features in our spectra could arise from the high- 
electron-concentration parts of the sample, where the en- 
ergy difference between adjacent Landau levels is small. 

Another factor possibly affecting the data could be 
lateral inhomogeneity within a single graphene plane, 
or even fragmentation of the layer - this could explain 
the weaker intensities observed for sample 2. Individual 
graphene planes in epitaxial graphite may be much more 
weakly coupled than is usually accepted for graphite. 
Several graphene layers 24| may, depending on the stack- 
ing scheme, exhibit linear and/or parabolic dispersion 
relations. Finally, we note that carriers with linear dis- 
persion may also be found at the H point of bulk graphite 
|25l I2H |27j . although the significant difference between 
the value of c found in [27J (0.91 xlO 6 m/s) and our ex- 
periment, as well as the structure of the samples, makes 
such an interpretation unlikely. The current experiment 
shows an absorption in good agreement with that ex- 
pected for graphene, but the simple approximation used 
does not exclude more complex scenarios. 

Concluding, we have measured the optical excitation 



spectrum of (relativistic-like) Dirac fermions in a con- 
densed matter system. These fermions are found in thin 
layers of epitaxial graphite, probably in single (or ex- 
tremely weakly coupled) graphene layers (or parts of lay- 
ers). Cyclotron resonance like transitions coexist with 
electron-hole (particle-antiparticle) like transitions, with 
energy positions and oscillator strengths in surprisingly 
good agreement with expectations based on a model of 
non-interacting particles with linear dispersion. 
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